Partial spin freezing in the quasi-two-dimensional La 2 (Cu,L«)04 
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In conventional spin glasses, the magnetic interaction is not strongly anisotropic and the entire spin 
system freezes at low temperature. In La2Cuo.94Lio.06O4, for which the in-plane exchange interaction 
dominates the interplane one, only a fraction of spins with antiferromagnetic correlations extending 
to neighboring planes become spin-glass. The remaining spins with only in-plane antiferromagnetic 
correlations remain spin-liquid at low temperature. Such a novel partial spin freezing out of a spin- 
liquid observed in this cold neutron scattering study is likely due to a delicate balance between 
disorder and quantum fluctuations in the quasi- two dimensional 5= 1/2 Heisenberg system. 



The parent compound for high transition-temperature 
superconductors, La2Cu04, is an antiferromagnetic (AF) 
insulator. Magnetic exchange interaction J between 
the nearest neighbor 5=1/2 spins of Cu 2+ ions in the 
Cu02 plane is several orders of magnitude stronger than 
the interplane exchange interaction, making quantum 
spin fluctuations an essential ingredient for magnetic 
properties in the quasi-two-dimensional (2D) Heisenberg 
system[H|ll3. The Neel temperature T/v of La 2 Cu04 
is suppressed rapidly by 2—3% hole dopants such as Sr, 
Ba or Li^, 0, while it is suppressed with isovalent Cu 
substitution at a much higher concentration close to the 
site dilution percolating threshold0. The strong ef- 
fect of holes has been shown to be related to induced 
magnetic vortices, which are topological defects in 2D 
systems 0- The paramagnetic phase exposed by hole 
doping at T <C J /kg is dominated by zero-point quan- 
tum spin fluctuations and is referred to as a quantum spin 
liquid lj. The approach of the Neel transition toward 
zero temperature promises detailed predictions for spin 
dynamics by extending the celebrated theory of critical 
phenomena for classical statistical systems to quantum 
statistical systems yj,|S|. 

However, in a wide doping range of La 2 Cui_ 2: Li a; 04 
below ~10 K, a spin-glass transition has been reported 
in muon spin rotation (/^SR)[f|, nuclear quadrupole res- 
onance (NQR);9] and magnetization^ studies. A sim- 
ilar magnetic phase diagram has also been reported for 
La 2 _ a; Sr a; Cu04 and Yi-xCaxBaaCusOeHl EI E3- In 
conventional spin glasses, magnetic interactions are more 
or less isotropic in space, and the whole spin system is 
believed to be frozen in the spin-glass phase[14j. If the 
spin-glass phase in hole-doped cuprates behaved as in 
conventional spin-glasses, the ground state would be a 
spin-glass instead of Neel order or quantum spin liquid, 
and as pointed out by Hasselmann et al.|15|. the quantum 
critical point of the AF phase would be preempted. 

Recently, 2D spin fluctuations in La2Cui_ a; Li a ;04 
(0.04 < x < 0.1) were observed to remain liquid-like 



below the spin-glass transition temperature Tgjl^, Il7| . 
The characteristic energy of 2D spin fluctuations satu- 
rates at a finite value below ^50 K as expected for a 
quantum spin liquidQ, instead of becoming zero at T g 
as for spin-glass materialsj0,^3- To reconcile these ap- 
parently contradicting experimental results, we have con- 
ducted a thorough magnetic neutron scattering investi- 
gation of La2Cuo.94Lio.06O4 to search for spin-glass be- 
havior. We found that in addition to the liquid-like 2D 
dynamic spin correlations, the rest of spins which par- 
ticipate in almost 3D and quasi-3D correlations become 
frozen in the spin-glass transition. This partial spin freez- 
ing in the laminar cuprate is distinctly different from to- 
tal spin freezing in conventional 3D spin-glass materials. 
The observed phase separation into spin glass and spin 
liquid components of different dimensionality sheds light 
on a long-standing confusion surrounding the magnetic 
ground state in hole-doped cuprates. 

The same sin gle crystal sample of La2Cuo.94Lio.06O4 
used previously [16j was investigated in this work. 
T g «8 K was determined using /iSR 5]. The lattice 
parameters of the orthorhombic Cmca unit cell are 
a=5.332A, 6=13.121 and c=5.402A at 15 K. Wave-vector 
transfers q near (000) and (100) in both the (hOl) and 
(hkO) reciprocal planes were investigated at NIST using 
the 30 meter high resolution small angle neutron scatter- 
ing (SANS) instrument at NG7, and cold neutron triple- 
axis spectrometer SPINS. We set the array detector of 
NG7-SANS to 1 and 9 m, corresponding to a q range 
from 0.012 to 0.39 A" 1 and from 0.0033 to 0.050 A -1 , 
respectively. At SPINS, the (002) reflection of pyrolytic 
graphite was used for both the monochromator and an- 
alyzer. Horizontal Soller slits of 80' were placed before 
and after the sample. A cold BeO or Be filter was put 
before the analyzer to eliminate higher order neutron in 
the fixed Ef=3.7 or 5 meV configuration, respectively. 
Sample temperature was controlled by a pumped 4 He 
cryostat which could reach down to 1.5 K. 

Hole induced ferromagnetic exchange has been theo- 



2 



y 1 1 




□ 3K ' 




J 


1 


W 


H 








X Ml 








i . 


.. i .... i 



0.01 0.02 03 

q (A"') 

FIG. 1: Measured SANS cross-section in the 10 pixels wide 
rectangle shown in the inset at 3 K (squares) and 30 K (dia- 
monds), and the difference of the intensity between 3 K and 
30 K (circles) as a function of wave vector transfer q. Lines 
are guide to the eye. Inset: the SANS pattern at 30 and 3 K. 
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FIG. 2: Magnetic quasielastic and elastic scattering along a) 
an in-plane direction and b) the interlayer direction at various 
temperatures. Open squares in b) was measured at (1.06,fc,0) 
and represent background. The solid lines are resolution con- 
voluted S 3D {q, E) + S q3D (q, E) in Eq. 



retically proposed in the Cu0 2 planeHHli. Although 
long-range ferromagnetic order has never been observed, 
there is the possibility of short-range ferromagnetic spin 
clusters which freeze in the spin-glass state. A powerful 
tool to probe such clusters is SANS|18|. Two reciprocal 
zones were studied, with incident beam parallel to the 
(001) or (010) direction, to cover all spin orientations. 
The experiments were carried out at 3, 10, 15, 30 and 
80 K. A collection time of 1 or 2 hours per temperature 
provides good statistics. No temperature dependence in 
the scattering patterns could be detected. The inset to 
Fig. ^ shows SANS patterns at 3 K and 30 K with inci- 
dent beam parallel to the (001) direction. Intensity at 3 
and 30 K in the rectangular box on the SANS pattern is 
shown in the main frame. The difference intensity (cir- 
cles) fluctuates around zero, and its standard deviation 
sets a upper limit of 1.5xl0~ 7 bn or 1.4xl0 _3 /zs per Cu 
for ferromagnetic moments in the clusters. 

While no appreciable ferromagnetic signal was de- 
tected for La2Cuo.94Lio.06O4, antiferromagnetic scatter- 
ing was observed along the rods perpendicular to the 
CuC>2 plane and intercepting the plane at the com- 
mensurate (7T, 7r)-type Bragg points of the square lat- 
tice. This means that antiferromagnetic correlations 
in the CuC>2 plane is chessboard-like, which is sim- 
pler than in L^-^Sr-cCuC^ 20]. Elastic and quasielas- 
tic scans through such a rod at (100) are shown in 
Fig. Ufa) at various temperatures. Inelastic scans have 
been presented in ^(|. There is little change in the 
peak width in these scans, consistent with previous 
results of temperature independent in-plane correla- 
tion length for La2Cuo.95Lio.05O4 2] and La2- a; Sr : j;Cu04 
(0.02< x <0.04)[2l| below 300 K. Modeling the width of 
the rod in La2Cuo.94Lio.06O4 with Lorentian 



^ > 274A from scan at (100) and E = 0, > 70l 
at (1,0.4,0) and E = 0, and ^ D > 55A at finite E from 
deconvolution. The □ denotes the in-plane component. 
These planar AF clusters correlate in three different ways 
in the interlayer direction, giving rise to almost 3D, quasi- 
3D and 2D magnetic correlations. Let us now examine 
the three components. 

Fig- 13b) shows elastic and quasielastic scans along the 
rod in the interlayer direction at 1.5, 49 and 180 K. Mag- 
netic intensity is composed of broad and sharp peaks at 
magnetic Bragg points (100) and (120) of the parent com- 
pound. The (110) peak is temperature-independent thus 
nonmagnetic. Fitting the broad peaks to Eq. JIJ, we ob- 
tained an interlayer correlation length £ q3D = 6.2(4) A. 
Again, no temperature dependence can be detected for 
£q3D De i ow 49, k Thus, the quasi-3D spin correlations 
are typically three planes thick. For the sharp peak, 
£ 3D > 168A Therefore, the number of correlated AF 
planes is more than 50, resembling 3D AF order. Both 
the broad and sharp peaks are energy-resolution-limited 
with the half-width-at-half-maximum < 0.07 meV, as it 
is displayed in Fig. [3] and in [l^. However, they should 
not be regarded automatically as from static magnetic or- 
der. Static magnetic signal was observed in AiSR study 5] 
only below T g =8 K at the spin glass transition. Thus, the 
energy spectra for the quasi-3D and almost 3D correla- 
tions may be modeled by £ 1 / e (i?), with e/h much smaller 
than 17 GHz=0.07 meV/h, the frequency resolution at 
SPINS, and larger than the zero- field /iSR static cutoff 
frequency of about 1 MHz H 03 for T > 8 K. 

The 2D AF correlations have been investigated in de- 
tail in [16j. The dynamic magnetic structure factor, 



W + (q0 2 ], (1) 



has been measured for E < 4.2 meV between 1.5 and 



FIG. 3: Measured S(q, E) as a function of E and interlayer 
k at 1.5 K, showing three color-coded magnetic components 
in Eq. (5). The E — peak is energy-resolution limited. 



FIG. 4: Temperature dependence of spectral weights I 
(green), I q3D (blue) and I 3D (red). See text for details. 



150 K near a magnetic Bragg wave- vector r, and k = 
q t. The almost 3D and quasi-3D spin correlations 
can be described, respectively, by 

5 3D (q, E) = I 3D ]T £$? (n a )C^ (k - r k )jO}/ e (E) (3) 



and 

S q3D (q, E) = I q3D <&°° (*o)£ $q3D (* - T k )C 1/c {E). 

T 

(4) 

With negligible ferromagnetic correlations, the total dy- 
namic structure factor is a summation of Eq. ©-0, 

S(q, E) = S 2D (q, E) + S q3D (q, E) + S 3D (q, E). (5) 

Of the four variables of 5(q, E), q^ is fixed at the (tt, ir)- 
type Bragg points. To comprehend the composition of 
S(q,E), it is sufficient to plot S(q,E) as a function 
of E and the interlayer wavenumber k. Such a plot of 
measured S(q,E) at 1.5 K is shown in Fig. EI with the 
temperature and q independent incoherent scattering at 
E = subtracted. The sharp peak fitted by the red 
curve is from S 3D (q,E), the narrow blue ridge at E=Q 
from S q3D (q, E), and the green surface from S 2D (q, E). 
The red peak at (100) is about one/three orders of mag- 
nitude stronger than the peak intensity of the blue/green 
surface, respectively. 

The spectral weights / dqdE S 3D (q, E)=I 3D and 
/ dqdE S q3D (q, E)=I q3D can be obtained by fitting 
resolution-convoluted Eq. ©-0 to scans in Fig. ©b). 
They are shown as a function of temperature in Fig. 0] 
with I 3D magnified by a factor of 5 for clarity. For the 
2D component, 
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Green circles in Fig. 0] represent I 2D with energy inte- 
gration from —0.5 to 4.2 meV, the same interval as in 



experiment 16] . Green squares are obtained with a wider 
integration range, \E\ < 10 meV, using the analytical ex- 
pression of x"{E) m El to extrapolate to E=10 meV, 
where spin fluctuations were observed in La2Cuo.9Lio.1O4 
using a thermal neutron spectrometer[22j. Thus, the 
green symbols serve as the lower bound of the 2D spectral 
weight which has ±00 as the integration limits. 

I 3D and I q3D appear simultaneously below ^150 K. 
Their concave shape in Fig. 01 differ drastically from the 
usual convex-shape of a squared order parameter, which 
was observed in ^SR study below T g =8 K 5]. They 
are typical neutron scattering signal from slow dynamic 
spin correlations in spin-glasses [3 which fluctuate 
in the frequency window between 1 MHz and 17 GHz 
for T > 8 K, and below 1 MHz for T < 8 K. Previ- 
ously, energy-resolution-limited neutron scattering from 
Lai.94Sro.o6Cu04 was observed to have a similar temper- 
ature dependence as I 3D in Fig. 01 and was attributed to 
spin freezing^] • The kink at 20 K of I 3D reflects an 
increased T g from 8 K to 20 K when_probing frequency 
is increased from 1 MHz to 17 GHz[ m ll4 | . 

The fact that I 3D decreases below T ff (17MHz)«20 K 
while I q3D continues to increase indicates that the 
"Edwards- Anderson order parameter" 0, |2^| dis- 
tributes only along the (1/cO) line. In conventional 
spin-glasses, the "Edwards- Anderson order parameter" is 
more isotropically distributed in the a-space [l4l Hal l23j . 
Thus, the spin-glass state in La2Cuo.94Lio.06O4 is charac- 
terized by interlayer disorder which upsets phase correla- 
tion between large AF clusters in different Cu02 planes. 
Another important difference from conventional spin- 
glasses in which all spins freeze below T g is that only a 
fraction of spins freeze in La2Cuo.94Lio.06O4. Other spins 
in 2D correlations remain fluctuating down to 1.5 K. This 
is consistent with previous theory that quantum fluctua- 
tions preve nt sp in-glass transition for 2D 5=1/2 Heisen- 
berg system 24]. The spin-glass component in our sample 
has to acquire interlayer correlations to achieve a higher 
dimension in order to be realized. It appears that the 
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lower critical dimension for a S—l/2 Heisenberg quantum 
spin glass is between 2 and 3. A further difference from 
conventional spin-glasses, for which one can measure the 
narrowing of magnetic spectrum toward £^=0|l8|. is that 
when S 3D {q,E) and S q3D (q,E) in La2Cuo.94Lio.06O4 
become detectable at about 150 K, they are already 
energy-resolution-limited. This property of S 3D (q, E) 
and S q3D (q, E) resembles the classic central peak phe- 
nomenon in the soft phonon transition|25l|. The dis- 
parate dynamics of the central peak and phonon are ex- 
plained by Halperin and Varma|2(| using a phase sepa- 
ration model: defect cells contribute to the slow relaxing 
central peak while coherent lattice motions (phonons) to 
the resolved inelastic channel. This mechanism has been 
applied with success to a wide class of disordered relaxor 
ferroelectrics|27ll28|. 

For La2Cuo.94Lio.06O4, we envision that disorder ac- 
companying doping prevents the long-range order of the 
AF phase mainly by upsetting interlayer magnetic phase 
coherence, see Fig.|3for the q-distribution of frozen spins. 
This upsetting is not uniform in the Griffiths fashion [2^ 
with weak and strong coupling parts in the sample. In 
our laminar material, however, the weak and strong cou- 
pling parts have different dimensionality: 2D and nearly 
3D, respectively. The 2D part is a spin liquid and repre- 
sents essentially the whole system at high temperature, 
see Fig. 0] Part of sample with stronger intcrplane cou- 
pling tends to order three dimensionally below ^150 K, 
producing S 3D (q,E) and S q3D (q, E) . The condensation 
of the 2D spin liquid at ~ 150 K into the 3D dynamic clus- 
ters of diminishing energy scale, instead of a true long- 
range order, may reflect the divergent fluctuations which 
destabilize static order below the critical temperature for 
random XY or Heisenberg systems |30J. The nearly 3D 
spin-glass instead of a 3D antiferromagnet finally orders 
at a much reduced T g » 20 K, when I q3D + I 3D ap- 
proaches the 2D spectral weight (Fig.0J. The phase sep- 
aration into spin liquid and spin glass components may 
be a consequence of no true mobility edge separating fi- 
nite and infinite range correlations at 2D|3J|. It would be 
interesting to extend the nonrandom theories^] to incor- 
porate this disorder scenario and to explain the observed 
crossover from the E/T to a constant E scaling in the 
2D spin liquid[TE[l3. 

In summary, spins in La2Cuo.94Lio.06O4 develop dy- 
namic AF order in the CUO2 plane with very long fp 
below 150 K. The characteristic energy of the 2D spin 
fluctuations is 0.18fc B T for T > 50 K and 1 meV for 
T < 50 K0. Below ^150 K, interlayer phase coherence 
appears between some of these planar AF clusters with 
an energy scale smaller than 70 /xeV. While the 2D AF 
correlations in an individual plane remain liquid down to 
1.5 K, coherent multiplane AF correlations become glassy 
below T g . The phase separation into 2D spin- liquid and 



"3D" spin-glass is most likely related to quasi-2D nature 
of magnetic exchange in the cuprates and is distinctly 
different from conventional spin-glasses. The theory of 
the spin-glass phase in doped cuprates should consider 
interlayer coupling, despite it is weak. Investigation of 
disorder effect on quantum critical spin fluctuations is 
called for. The phase separation, instead of a uniform 
magnetic phase, suggests a possibility that superconduc- 
tivity and the almost 3D AF order reside in different 
phases in La2- 2: Sr a; Cu04 and Yi_ a Ca a Ba2Cu306+j / . 
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